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Abstract. We propose to study echo dynamics in a random matrix framework, 
where we assume that the perturbation is time independent, random and orthogonally 
invariant. This allows to use a basis in which the unperturbed Hamiltonian is diagonal 
and its properties are thus largely determined by its spectral statistics. We concentrate 
on the effect of spectral correlations usually associated to chaos and disregard secular 
variations in spectral density. We obtain analytic results for the fidelity decay in the 
linear response regime. To extend the domain of validity, we heuristically exponentiate 
the linear response result. The resulting expressions, exact in the perturbative limit, 
are accurate approximations in the transition region between the "Fermi golden rule" 
and the perturbative regimes, as examplarily verified for a deterministic chaotic system. 
To sense the effect of spectral stiffness, we apply our model also to the extreme cases 
of random spectra and equidistant spectra. In our analytical approximations as well as 
in extensive Monte Carlo calculations, we find that fidelity decay is fastest for random 
spectra and slowest for equidistant ones, while the classical ensembles lie in between. 
We conclude that spectral stiffness systematically enhances fidelity. 



PACS numbers: 05.45.Mt, 03.65.Yz 



Fidelity decay 



2 



1. Introduction 

Loschmidt echoes, proposed more than a century ago as a gedanken experiment P , have 
been realized experimentally j2] measuring the corresponding auto correlation function 
known as fidelity Fidelity has also been used as the simplest benchmark for the 
reliability of quantum information devices |1] , and recently Kaplan ^ has used fidelity of 
eigenstates to test uncertainty of quantisation. Many papers study general properties of 
fidelity under different regimes, including the semiclassical jHl Ej , and the linear response 
regime j^lHUHl- The latter can be divided into the so called (standard) perturbative and 
the "Fermi golden rule" regimes ^Oj- In the perturbative regime, linear response theory 
will be correct for times long as compared to the Heisenberg time, while in the golden 
rule regime, it will break down before the Heisenberg time. The basic aim of this paper 
is to use a random perturbation V and, following the spirit of random matrix theory, 
to derive the universal features we may expect to be relevant for chaotic systems. The 
randomness of the perturbation implies orthogonal (or in the time reversal breaking 
case unitary) invariance. This in turn permits to use a basis in which the unperturbed 
Hamiltonian Hq is diagonal. Then, only the spectral properties of Hq enter the problem. 
We disregard secular variations in the spectral density and concentrate on fluctuations. 

Our treatment shall make extensive use of linear response theory, because this 
will cover the basic needs of quantum information. Yet the model is not restricted to 
the lowest orders in this framework, and we shall see in Monte Carlo simulations how 
to extend heuristically our description to long time scales, exponentiating the linear 
response result. This will prove exact in the perturbative limit and will give good 
results even for moderately small perturbations. Note that in any case, we expect the 
linear response theory to hold much longer than the time corresponding to the inverse 
spectral span for which the eigenphases of the time evolution operator fill a small part of 
the unit circle, and which we call Zeno time. This regime allows for the straightforward 
expansion of the exponentials and thus is of limited interest. Furthermore in the context 
of information processing the linear response regime is most relevant, because in the 
semiclassical regime it will be almost impossible to have coherence in the quantum 
registers, and times short compared to the inverse spectral span are typically too short 
for quantum information processes. 

We therefore construct the following model: Hq is a diagonal matrix with an 
unfolded spectrum (constant level density) and is a random matrix pertaining to the 
Gaussian orthogonal or unitary ensemble (GOE or GUE). To model chaotic systems we 
shall follow the quantum chaos conjecture jlTj, and use unfolded GOE or GUE spectra 
for Hq. To contrast these with the possible extreme cases, we shall also consider picket- 
fence (equally spaced) and random spectra. Actually, any spectrum can be used, whose 
form factor is known. Random spectra are often associated with integrable systems |12j . 
but we must keep in mind that the effects we describe in echo-dynamics are state 
dependent, and thus other important properties of integrable systems such as invariant 
tori will typically also be relevant. We can therefore expect the GOE and the GUE 
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cases to relate directly to chaotic systems, and indeed we shall see that such a relation 
can be established, as shown in ^21- Indeed, the random matrix model used there, is 
very similar to the one we propose. Note that the use of a random perturbation has 
been proposed in reference |3]. 

The other important aspect in echo dynamics is the dependence on the initial 
state. The analysis of coherent states - so important in other studies - is not considered 
here, because, typically, we do not expect to have an underlying classical model at our 
disposal. The most important property is then the spectral span of the initial state. We 
shall mainly consider two extremes: the evolution of an eigenstate of the unperturbed 
system, and the evolution of a random state with a given, but fairly large spectral span. 
The latter is the relevant case for quantum computing. This is easily understood if we 
think e.g. of Fourier transforms. The transforms of simple functions such as Gaussians 
are readily done analytically, while we would like to exploit the quantum algorithm for an 
arbitrarily complicated case. Yet it is also interesting to consider eigenfunctions partially 
because of the applications in the context of effects of different types of quantisation [H] 
but also because these display most markedly the effects of spectral fluctuations. 

We measure the distortion of a state by calculating the fidelity, which may be 
interpreted as the autocorrelation function of a forward evolution with Hq and a time 
reversed evolution with the perturbed Hamiltonian H = Hq + XV. The same quantity 
can also be interpreted as a cross correlation function of the same initial state evolving 
under the two Hamiltonians Hq and H. Specifically, if and 'ifo{t) are the functions 
evolving under H and Hq, respectively, from the same initial function \l/(0) = \l/o(0), we 
define the fidelity amplitude as the matrix element 

m = {Mt)\m) = imiuoi-t) u{t)\m) , a) 

where Uo(t) is the unitary propagator associated with Hq, and U(t) is the one associated 
with H. From the point of view of a quantum Loschmidt echo, UQ{—t)U{t) would 
be the "echo"-operator. The fidelity is usually defined as F(t) = In the 

interaction picture with the state of the system denoted by x(t), the fidelity is simply 
the autocorrelation function. To see this we recall that 

^(t) = Uo{t) x{t) ihdt x{t) = XV{t) x{t) , (2) 

where V{t) = Uo{-t) V Uo{t). Hence: 

fxt) = {xmuo{-m{t)) = {xmm ■ (3) 

In the next section we recall the relevant linear response relations in the interaction 
picture, and apply them to the calculation of the fidelity amplitude. In the perturbative 
regime, i.e. for small A, the Born series can be summed to yield the well known 
Gaussian decay. We therefore give a phenomenological formula for the general case 
by simply exponentiating the linear response result, and we compare to Monte Carlo 
simulations to establish the accuracy of this formula. In section |23l we apply our results 
to the standard map, using calculations by Cerruti and Tomsovic |T2]. We find that our 
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heuristic formula is not only simpler, but also closer to the numerical experiment than 
the expression derived there. 

In section 01 we proceed to calculate the fidelity. We observe a marked difference 
between the behaviour of a single eigenstate and that of a superposition. In the former 
case, the differences resulting from different spectral statistics are enhanced. Again we 
compare with numerical calculations. In section |3] we study the asymptotic behaviour 
of the fidelity for long times in various situations, and we obtain some analytical results 
from the perturbation theory. All these studies use spectral correlation functions in the 
large limit for the GOE (and GUE) case. Finally we give some concluding remarks 
in sectional 

2. Fidelity amplitude 

In this section we shall first derive a linear response formula for the decay of the 
fidelity amplitude. The central quantity of interest is here the correlation integral, 
which is a double integral over the autocorrelation function of the perturbation in 
the interaction picture. Given that the perturbation is taken from the GOE (GUE), 
the correlation integral only depends on the spectral properties of Hq, and we shall 
evaluate the correlation integral for different cases. Note that, unless stated otherwise, 
GOE perturbations are used throughout this article. We then discuss the possibility to 
extend the validity of the linear response result by packing it into an exponential. This 
phenomenological formula works very well up to the crossover between the perturbative 
and the Fermi golden rule regime. In the last part of this section we discuss the sample 
fluctuations of the fldelity amplitude concentrating on the regime of high fldelity, as this 
is the one most relevant for quantum information applications. 

As a preliminary step, let us introduce more convenient units for time and energy. 
In order to take advantage of the orthogonal (unitary) invariance of the perturbation, 
we write the Schrodinger equation in the eigenbasis of Ho: 

Ho = diag(^„) , f/o(t) = diag(e-^^"*/^) . (4) 
Thus, we obtain in the interaction picture: 

ih dt x{t)=X diag(e'^" V diag(e-'^^ x{t) . (5) 

Denoting with d the average level spacing in the spectrum of Ho (for simplicity we shall 
assume that d is constant in the energy range of interest), the Heisenberg time is deflned 
as tn = 2'Kh/d. It is convenient to measure time in units of tn and energy in units of 

d, introducing the dimensionless time s = t/tu and spectrum {e^ = E^/d}. Then we 
obtain for x'{s) = x{tH s): 

/.fr A ° ° 

ids x\s) = — - diag(e2"'""^) V dmg{e-^^'^^') x\s) . (6) 

This shows, that without restriction of generality, we may assume that the Heisenberg 
time and the average level spacing are both equal to one - and in the rest of the paper 
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we shall indeed do so. This allows us to stick to the symbols for time and energy, as 
introduced in the beginning. 

2.1. Linear response 

o _ 

With Uo{t) = diag[exp(— 27ri Eq, t)], and the shorthand V{t) = Uo{—t)VUo{t), we can 
approximate x{t) using the Born series: 

=x(0) -27riA /" dr V(r) x("-^)(r) , = x(0) . (7) 



JO 

Up to second order we obtain: x^^^t) = X{t) x(0), where 

X{t) = l-2mX dr 1/(r) - 47r2 AM dr / dr' V{t)V{t') . (8) 
Jo Jo Jo 

Note that X{t) is the linear response approximation of the echo operator, mentioned 

above. In general, we are interested in the various moments of the components of x^'^^{t). 

Their knowledge will allow us to calculate averages of the fidelity, as well as the survival 

probability or even the purity. However, we may additionally choose the initial state 

x(0) to be random and eventually even Hq could be a member of some ensemble. In 

particular, if we want to keep the initial states arbitrary, it is more convenient to consider 

the various moments of the matrix elements Xapit), when the average is performed over 

the GOE (or GUE) matrix V. We define these matrix ensembles such that 

/T. T. \ / ^^l^^k • GUE 

[ OikOji + Oiidjk : GUE 
Our linear response approach is very similar to the one by Prosen and coworkers |14j . 
However, in their work, Hq and V are fixed while the average is taken over initial 
conditions only. Here we will primarily average over the perturbation V, thereafter over 
Hfj, if we consider an ensemble, and finally over initial states, if these are random. 

The spectral correlation function In order to average the linear response operator X(t), 
given in (jH)), we basically need the correlation function: {V{r)V{T'))v, where {■ ■ .)v 
denotes the ensemble average over the random matrix V. Here, we present the full 

o 

calculation for the GOE perturbation, only. With Aa{t) = exp(— 27ri E'q, t), we obtain: 
{[Vir) ViT')U)y = 5^ A„(-r) (K, A,(r - r') V,,)y A^(r') 

7 

= ^c.{-t) A^(r - t') A^(r') [S^^S^p + Sa/s] 

7 

= 5a/3C,(r'-r), C7,(t) = l + 5^e2-'(^--^")*. (10) 

7 

The correlation function can be expressed in terms of the two-point form factor 
b2(t), one of the important fluctuation measures in quantum chaos studies jT^j. Namely, 
if we denote with (. . .)o the average over the diagonal matrix Hq, we may write: 

{C^{t))o = 2 + 6{t)-b2{t) . (11) 
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The i^o-average makes sure to obtain a well behaved, i.e. smooth, two-point form 
factor. Though, note that in practice, this additional averaging may be avoidable. The 
corresponding expression for the GUE perturbation is: 1 + 6{t) — 62 (^)- 



Let us now turn to the calculation of the average echo operator in the linear response 
approximation: 

(X„^(t))y = - Atx^ X' [ dr [ dr' C„(r') . (12) 

JO JQ 

Additional averaging over Hq gives a scalar result: 

((^a/3(t))y)o = l-47r2A2C(t), (13) 

where 

C{t) =e + t/2 - f dr r dr' 62 (r') . (14) 
Jo Jo 

Note that the integration of the delta function gives |. For the GUE perturbation, the 
t^-term should be cut in half. Due to the scalar result in equation |T!?jl . the average 
fidehty amplitude does not depend on the initial state: 

{{m)v)o = {{{X{t))y)o)o = 1 - 47r2 C{t) . (15) 

Here (. . .)o is a short hand for the expectation value with respect to the initial state. If 
we avoid averaging over Hq, the average fidelity amplitude reads: 

{f{t))y = {{X{t)).)y = l-4n'X' fdr F dr' V|x,pC,(r'). (16) 

-'0 Jo „ 

If the participation ratio of the initial state x is sufficiently large, i.e. if the inverse 
participation ratio, ipr(x) = Yla I 

X Q I ^ IS sufficiently small, this may lead to self- 
averaging - just as it does for the autocorrelation function of certain deterministic 
systems. Then we may still use the general result, equation ()15p. inserting the proper 
two-point form factor for this case. 

2.2. Correlation integral 

Clearly, the correlation integral C{t) is the essential quantity which determines the 
fidelity (amplitude) decay. It is shown in figure [T] for the Poisson, GOE, GUE, and 
picket-fence spectrum. The two-point form factor b2{t) is zero in the Poisson case, while 
for all other cases, it is given in Appendix A The level repulsion, present in all but 



the first spectrum, tends to slow down the increase of the correlation integral. This is 
simply a consequence of the so called "correlation hole" ^H] , which can be found in the 
integrand of equation (fT^ . Note however that for large times the quadratic increase 
with time remains essentially unaffected [panel (a)]. In figure ^b), we subtract the 
t^-term in order to display more clearly the particular effects of the different spectral 
correlations. While in the Poisson case the remaining term is linear, in the correlated 
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Figure 1. (a) The correlation integral C{t), equation H14|) . for the Poisson spectrum 
(solid line), the GOE spectrum (dashed line), and the picket-fence spectrum (dotted 
line), (b) The function C{t) — for the three spectra in panel (a), and for the GUE 
spectrum (dash-dotted line). 



spectra (GOE, GUE and picket-fence) the linear increase at the origin is capped. In 
fact, we can calculate explicitely the asymptotic behaviour of the correlation integral. 
In the GOE case, using the exact formula ()A.11|) . we obtain: 

C{t) =t +2~ J,^ = + 12 + ^^^"^ • ^ 

In the picket-fence case, the correlation integral (with the term subtracted) has zeros 
at integer values of t. The GUE spectrum leads to C{t) — approaching 27r^/3 from 
below. 



2.3. Exponentiated linear response 

In the perturbative regime, i.e. in the limit of small A, the fidelity decay is known to be 
Gaussian [3]. This can readily be reproduced in the present context. The Born series 
can be trivially summed under the assumption that eigenstates are not affected by the 
perturbation. The result corresponds to the exponentiation of equation ()15|) : 

(/(t))=exp[-4vr2A^C(t)]. (18) 

This result is exact in the limit 27rA ^ 1, only, but we shall show below that errors 
are fairly small even for A ~ 0.1, and negligible for A ~ 0.01. The essential point is 
that the correlation integral, whether used in equation (jl5|) or (jl8p clearly displays the 
transition between a linear decay at short times and a quadratic decay at long times. 
Accordingly, for small A, the decay of the fidelity amplitude is dominantly Gaussian, 
whereas for larger A, it is dominantly exponential [HE]- 

Figure |21 shows the behaviour of the fidelity amplitude according to equation (fTSj) 
for the GOE and the Poisson case, in the large N limit. For the latter, the linear and 
quadratic terms are also shown separately. Note that the GOE results lie systematically 
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Figure 2. The fidelity amplitude in the crossover regime, from linear to quadratic 
decay. 27rA = 1. Thick solid and dashed lines: Theory for the Poisson and the 
GOE case respectively, i.e. equation 1)15(1 packed into an exponential (see text for 
details). Thin dotted lines: Theory for the Poisson case, where either the linear term 
or the quadratic term has been ignored. The inset shows a zoom into the area around 

t = 0, (/W) = i. 
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Figure 3. Difference between the Monte Carlo and the theoretical result, 
equation (|18|l . for the fidelity amplitude, for the Poisson (solid line), the GOE (dashed 
line) and the picket-fence (dotted line) case. In (a) A ~ 0.1, whereas in (b) A — 0.01 
(note the different scale on the ordinate). 



above the Poisson ones due to the correlations in the spectrum of Hq. The crossover 
region is precisely the one of most interest and, therefore, we check the accuracy of the 
phenomenological expression (fTH|) beyond the perturbative limit. 

Figure El shows the deviation of the fidelity amplitude obtained for Monte Carlo 
calculations for = 100 from what we find using equation pSj) and the = oo 
correlation integral. Besides the GOE case, we considered also spectra without 
correlations (Poisson case) and a spectrum with equally spaced levels (picket-fence case). 
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In this way we cover a broad range of possible spectral correlations. For A = 0.1, 
figure a), the deviations are of the order 0.5% of the maximal value for the fidelity 
amplitude. We checked that these deviations are not due to finite size effects. By 
varying and also by considering the correlation integrals for finite N (for the Poisson 
and the picket-fence case), we found that those effects do not alter the error noticeably, 
except for a small "hump" ai t ^ 0.25 in panel (a). In panel (b), the error is shown 
for A = 0.01. Its absolute value is further reduced by a factor of ten, roughly, and 
finite size effects are no longer observable. Note that the errors for the random and the 
picket-fence spectrum are of the same order of magnitude as for the chaotic case. This 
suggests, that our exponentiated formula works well in all cases, where the two-point 
correlations are known. In order to compare the error curves in figure El with the actual 
value of the fidelity amplitude, note that the time thaif where the fidelity has dropped 
to half of its initial value is about thaif ~ 1-1 in panel (a), and thaif ~ 13 in panel (b). 
In all, we see that the results shown in figure |21 are very similar to what we can expect 
from a Monte Carlo calculation. 

2.4- Sample fluctuations 

In the figures Q and |2l we saw the effects of spectral correlations on the decay of the 
average fidelity. But are these effects observable also for an individual system? To 
answer this question, we shall now study the fluctuations of the fldelity amplitude, 
focusing on the regime of high fldelity. 

To obtain the numerical results shown in flgureEland hereafter, we average over the 
GOE perturbation V, the spectrum of Hq and eventually over initial states x. The initial 
state may be an eigenstate of Ho; then we choose it from the centre of the spectrum 
in order to minimise border effects. Alternatively, it may be a random (orthogonally 
invariant) state. Though the average value of the fldelity amplitude does not depend on 
the choice of the initial state, higher moments (such as the variance) do. 

Figure ID shows the average fldelity amplitude (/(t)) in the perturbative regime 
for the Poisson and the GOE case. We performed 10 independent ensemble averages, 
each of them over rirun = 1000 random systems. From this we calculated the statistical 
uncertainty for a single ensemble average. In this flgure, as well as in similar ones 
below, the area which is not more than one (two) standard deviations away from the 
total average over all samples, is plotted in dark (light) gray. In addition, we plotted 
the pure linear response result, equation (fT3j) for the Poisson ensemble (solid line) and 
the GOE (dashed line). While in panel (a), we show f{t) itself, in panel (b), the Poisson 
theory is subtracted in order to show more clearly the differences between the Poisson 
and the GOE case on the one hand, as well as between numerics and linear response 
theory on the other. 

The large sample-to-sample fluctuations are clearly due to the special choice of the 
initial state. In linear response theory, and after the averaging over V, it can be seen 
that the fldelity amplitude depends on only out of A^(A^ — 1) available eigenvalue 
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Figure 4. (a) Fidelity amplitude for N ~ 100, A = 0.01; the initial states are 
eigenstates of Hq. The thick solid and dashed lines show the theoretical result H15|l 
for the Poisson and the GOE case, respectively. The numerical results (shaded areas) 
are obtained from 10 small sample averages of size rirun = 1000, each. The total 
average of all samples lies in the centre of the gray bands, while their borders are 
given by plus/minus one standard deviation (dark gray), and plus/minus two standard 
deviations (light gray), (b) The same quantities as in (a), but with the Poisson theory 
subtracted. 
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Figure 5. Fidelity amplitude with the Poisson theory subtracted, for N = 100, A = 
0.01; random initial states. The numerical results (shaded areas) are obtained from 10 
small samples of size rirun = 1000, each. The significance of the gray scales is the same 
as in figure01 (a) The thick solid and dashed lines show the theoretical results for the 
Poisson case (here, we expect zero) and the GOE case, based on the linear response 
result p5(l . (b) The same data, but with the linear response result being packed into 
an exponential p8|l . 



differences [see equation |Tn|l ]. By consequence, we need very large samples in order to 
obtain accurate averages. However, the statistical ffuctuations can also be reduced by 
probing the decay of f{t) at the same system for various initial i^o-eigenstates, or by 
using a random initial state, as shown below. 

Figure El shows the fidelity amplitude (with the Poisson theory subtracted), for the 
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Figure 6. Fidelity amplitude for N 100, A = 0.01; small samples of size 
n-run = 1000. The mimerical results are represented by shaded bands, where the 
significance of the gray scales is the same as in figure 0] Data for the picket-fence 
spectra, where the initial states are eigenstates of Hq. The linear response result H15|l 
for the Poisson case (thick solid line) and the picket-fence case (thick dashed line). 

Poisson and the GOE case, where the initial states are taken to be random. We used the 
same sample size as in figure |3 Nevertheless, the statistical fluctuations have drastically 
decreased, and the Poisson- and the GOE-data can now be distinguished without 
difficulty. While in panel (a) we used the pure linear response theory, equation (fTH|l . 
for the Poisson and the GOE case, in panel (b), we used the exponentiated version, 
equation (jlSp . One can clearly see that there are systematic differences between the 
numerical data and pure linear response theory, while we obtained a near perfect 
agreement with the help of the exponentiation. 

Using a result from the following section El one can understand the reduction of 
the statistical ffuctuations. According to equation ()28|) . the variance of the fidelity 
amplitude is proportional to the IPR of the initial state. As for a random state 
X : ipr(x) = 3/(A^ -|- 2), it implies that the statistical deviations should be about 
^y (N + 2)/3 ~ 6 times smaller than in figure |3fb). We finally note that the error which 
remains after exponentiating equation is so small, that it will typically be lost in 
these fiuctuations. 

From figure El we can see that, also for the picket-fence spectra, the linear response 
theory does a good job. What is more interesting, is a comparison of the statistical 
deviations of the numerical curves for this figure, and figure l^a) which is practically 
the same plot, but for the Poisson and the GOE spectra. As the spectrum in the picket- 
fence case is no statistical quantity, it means that the statistical deviations are mainly 
due to the fiuctuations in the perturbation matrix V. This implies that the fidelity 
decay is sensitive to the level correlations in a single spectrum (of very moderate length, 
N = 100). One would just have to probe this particular system with different initial 
states, and different perturbations. Quantum dots with variable magnetic fields may be 
an appropriate test ground. 
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From the numerical analysis of the sample fluctuations, as well as from the linear 
response result for the variance of the fidelity amplitude, we conclude that the sample 
fluctuations are proportional to the inverse participation ratio of the initial state x. 
To observe a significant difference between a GOE spectrum and a spectrum without 
correlations, we need either very large samples (for initial i^o-eigenstates) , or initial 
states with very small inverse participation ratio {e.g. random initial states). The 
sample fluctuations are mainly due to the randomness of the perturbation V, not so 
much due to the random Hq. Finally, we have verified that the exponentiation of the 
linear response result improves the agreement with numerics considerably. In section 
to follow, we will perform a further, more realistic, test. 

2.5. A deterministic example 

In reference ^21 Cerruti and Tomsovic calculate the fidelity amplitude for the quantised 
standard map. They also present an analytical formula for the decay of the fidelity 
amplitude, which is based on a "uniform approximation". In essence, they use a 
random matrix model similar to ours, and an additional semiclassical argument which 
allows them to determine the perturbation strength. Their final formula agrees with 
ours in both limit cases of weak and strong perturbation, if we disregard the spectral 
correlations, i.e. if we set the two-point form factor in the correlation integral, 
equation ()14|). equal to zero. However, in the crossover region, there are further 
deviations. 

In order to compare our approach to theirs, we first have to adapt our model to 
a two-fold symmetry in the standard map, which leads to two independent symmetry 
sectors in the Hamiltonian. This infiuences the correlation integral, and by re-examining 
equation (fTUj) . we obtain: 



where the unit for time is again chosen such that = 1. Finally, we have to take 
into account that our strength parameter A, is related to Cerruti and Tomsovic's by: 
= 4A, and that in reference ^Hj the unit for time is chosen differently. In all, we 
obtain for the fidelity amplitude: 



where Cerruti and Tomsovic's semiclassical argument gives A ^ 0.0395. The same result 
follows also from general considerations of j3]. 

In figure U\ we reproduce the numerical data for the quantised standard map from 
Cerruti and Tomsovic, and compare it with a random matrix calculation and our theory, 
equation (|20|l . In the random matrix model we choose a matrix of dimension = 1000, 
which is block-diagonal. Each block contains a matrix of dimension A^/2, produced 
from an unfolded GOE spectrum (with average level spacing equal to two) and a 
GOE perturbation, as defined earlier. The variance of the non-zero off-diagonal matrix 




(19) 



fsymit) = exp[-16 7r^ A Csymit)] , 



(20) 
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Figure 7. (a) The fidelity amplitude for the quantised standard map (data from ll^p 
(solid line). Our exponentiated linear response result, equation (|20|1 for A = 0.0395 
(dashed line), and the theoretical result of Cerruti and Tomsovic (dotted line). The 
average over the full sample lies in the centre of the gray bands, while their borders 
are given by plus/minus one standard deviation (dark gray), and plus/minus two 
standard deviations (light gray), (b) The same curves as in (a) but with our theory, 
equation subtracted. 



elements is = 4A. To estimate the statistical error on the fidelity amplitude, we 
calculated f(t) for ten different random matrices of this type. In addition, we plot the 
corrected theoretical result from Cerruti and Tomsovic |13j . 

The fidelity decay for the quantised standard map is well described by our random 
matrix model. In particular in panel (b) we can see that the corresponding curve stays 
well within the statistical limits for individual members of the random matrix ensemble. 
Though our theory is not exact (cf. figure 01 and note that A ~ 0.4 is quite large), the 
systematic deviations from the random matrix results, as well as from the quantised 
standard map calculation, are very small. The theory of Cerruti and Tomsovic shows 
somewhat larger systematic deviations, in particular in the region around t = 0.2. Note 
that we compare with the semiclassical result of reference ^Hj, because we use the same 
classical parameters in our random matrix model. 

3. The fidelity 

In an experiment it is usually easier to measure the fidelity F{t) = \f{t)\^ rather than 
the fidelity amplitude. Though, often it is assumed that the average fidelity is simply 
given by the absolute value squared of the average fidelity amplitude, and we shall see 
under which conditions this is justified. In distinction to the average amplitude, the 
average fidelity depends on the choice of the initial state. For instance, if the initial 
state is an eigenstate of Hq, the fidelity F{t) coincides with the survival probability or 
autocorrelation function. 
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In order to calculate the fidelity -F(t), we first expand the linear response approximation 
of the echo operator (jH)) according to increasing powers of A: 

X{t) = 1 - 27ri A I{t) - A-n"^ A^ J{t) , 



m 



dr V{t) J{t) 



dr 



dr' V{t) V{t') . 



(21) 
(22) 



For an arbitrary initial state x, the average fidelity reads: 

{F{t))y = ( (X(t))o {X{t))l)y = l-4n'X' {2Re {{J{t)).)y - mt))^f)v) , (23) 

where we considered only terms up to second order. For the average ((</(^))°)\/ we 
obtain (cf. section ITTj) : 



{Jap{t))v = Sap fdr Tdr' C,(r - r') = S^p f dr T dr' C„(r') 
Jo Jo Jo Jo 



while 



{\{Iit))i?)v=Y.Vp°x,°x! dr [ dr' {Va^ 



p{r)V*{r'))y 



(24) 



(25) 



The average product of two matrix elements of V can be computed with the help of 
equation 0. For the GOE perturbation, we obtain: 

(K/3(r) V*{r'))y = [5a, A,(-r + r') A^(r - r') + A,(-r - r') A^(r + r')] 



o o 



o o 



5 5a p27ri(i?c-S/3)(r-r') , r e 2-Ki{Ec.-E p){t+t') 



^2TTi{Ec-Eft) (r-r') 



^27ri{Eo.-Ep) (t+t') 



ae "/37 

: a = P = ■y = e 

: a = 7 7^ /5 = £ 

: a = e ^ P = •y 
: otherwise 



(26) 



which leads to 



dr / dr' 



+5: 



° 1 2 I ° 



o o 



2 27ri(£;„-£;^) (t-t') 



o2 



o o 



. (27) 



Here = ipr(x) is the inverse participation ratio (IPR) of the initial state x. If 

the initial state is an eigenstate of Hq, then ipr(x) = 1, and this is the only term which 
survives. If x is complex, then we expect that term which depends on (r + r') to vanish. 
On the average, the contribution of each of the exponential terms is only of order N~^. 
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Figure 8. Fidelity for N = 100, A = 0.01; the initial states are eigenstates of 
Ho. (a) The thick solid and dashed lines show the linear response result (I28II for 
the Poisson case and the GOE case, respectively. The numerical results (shaded areas) 
are obtained from 10 small samples of size Tij-un = 1000, each. The picket-fence case 
is also included. The total average of all samples lies in the centre of the gray bands, 
while their borders are given by plus/minus one standard deviation (dark gray), and 
plus/minus two standard deviations (light gray), (b) As in (a) but for 10 large samples 
of = 50 000. In the picket-fence case, the theoretical result is plotted (thick 

dash-dotted line) instead of the numerical data. 



Therefore, the whole expression in the second hne of equation ()27|) is of order ^. In 
the hmit of large A^, it can be neglected, leading to 

{F{t))v = 1 - 2(27rA)2 |c(t) - ipr(x) t^j 

= |(/(t))vP + 2{2n\r ipr(x) + 0{t') . (28) 

For the GUE perturbation, equation (|26|) would simply read: {Vai3{T)V*^{T'))v = 

o o 

Sai3 exp[27ri(£'Q, — E13) (r — r')]. For the average fidelity, it would have the only 
consequence that all terms involving ipr(x) had to be multiplied by one half. 

The second line of equation (j^H|) gives the variance of the sample fluctuations of 
the fidelity amplitude, discussed in the previous section. It shows that, if the IPR of the 
initial state goes to zero, then indeed, the average fidelity is given by the absolute value 
squared of the fidelity amplitude (at least in the linear response regime). By contrast, 
if the IPR of the initial state is large (few principal components) , then the behaviour 
of the fidelity is quite different from that of the fidelity amplitude. In the extreme 
case that the IPR is equal to one (initial state is an eigenstate of Hq), the IPR-term in 
equation ()28p completely kills the t^-term in the correlation integral C{t). This holds 
equally true whether V is chosen from the GOE or the GUE. 

3.2. Average fidelity and sample-to -sample fluctuations 

In figure IHl we show the fidelity decay for states which are initially eigenstates of Hq. Due 
to the cancellation of the t^-terms, the fidelity decays much slower, and the effects of the 
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Figure 9. Fidelity for N = 100, A = 0.01; random initial states. Results for 10 
small samples of size ny-^n = 1000. The numerical results are represented by shaded 
bands, where the significance of the gray scales is the same as in figure |S1 (a) Data 
for the Poisson case, the GOE-case, and the picket-fence case. The thick lines give the 
theoretical results for the Poisson case (thick solid line), the GOE-case (thick dashed 
line), and the picket-fence case (dash-dotted line), (b) As in (a) but with the Poisson 
theory subtracted (linear response, packed into an exponential). 



spectral correlations are more pronounced. For example, for the picket-fence spectrum, 
we get practically complete revivals at integer multiples of the Heisenberg time, while in 
the Poisson case, the fidelity decays linearly. Note that exact revivals can occur in the 
perturbative limit, only. There the eigenvectors of H are arbitrarily well approximated 
by the eigenvectors of Hq. The statistical deviations seem to be much smaller in the 
picket-fence case; see figure ISI^a). However, this is simply because the fidelity is closer to 
one. We checked that the relative error for the correlation integral, remains comparable 
for all three cases. 

In figure IHfb) one can see a small difference between the theoretical linear response 
result and the numerics for both, the Poisson and the GOE case. This time, it does 
not help to exponentiate the linear response result. Here, the fidelity decay is too slow 
to produce an observable effect. Note also that the deviations differ in sign. While the 
theory underestimates the fidelity decay in the GOE case, it overestimates the decay 
in the Poisson case. This points at a more complicated behaviour of the fidelity (as 
compared to the fidelity amplitude). In the picket-fence case, by contrast, the difference 
between theory and numerical result is much too small to be resolved. Due to the 
fact that the statistical fluctuations have also diminished, we decided not to show the 
numerics at all. 

In figure El we show the behaviour of the fidelity for random initial states. As 
expected, the fidelity behaves very similar to the fidelity amplitude in this case {of. 
figures m and El) • Its decay is again dominated by the t^-term, and the effect of the 
correlations in the spectrum of Hq has decreased, as compared to figure |H1 Again, it 
makes sense to exponentiate the linear response result, in analogy to equation ()18|1 . 
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Figure 10. Fidelity for N = 100, A = 0.1; random initial states. Results for 10 
small samples of size n^un = 1000. The numerical results are represented by shaded 
bands, as in figure |S1 (a) Numerical results for the Poisson case, the GOE-case, and 
the picket-fence case. The thick solid line shows the pure linear response result for the 
Poisson case, (b) as in (a), but with the Poisson theory subtracted (linear response, 
packed into an exponential). The theoretical results, using equation (|29() . are shown for 
the Poisson case (thick solid line), GOE-case (thick dashed line), and the picket-fence 
case (thick dash-dotted line). 



Figure ITUl shows the fidehty decay for random initial states in the crossover regime, 
A = 0.1, between perturbative and golden rule decay. According to the discussion of the 
figures d and |21 we expect a stronger effect of the spectral correlations in this regime. 
Note e.g. at the Heisenberg time, t = 1, the relative differences between Poisson and 
GOE as well as between GOE and picket-fence, are of the order of 10%, each. The 
theoretical curve in figure [TUT a) is the pure linear response result. It is obviously not 
sufficiently accurate. To improve, we use exponentiation, but in addition we also take 
into account a rather trivial finite size effect: Namely, as will be shown in the following 
section: limt-too F{t) = Coo ~ 3/(A^-l-2). Therefore, instead of the simple exponentiation, 
we use 

h{x) = (1 - Coo) e-^/^i-'^-) + Coo , Coo ^ ■ (29) 

This phenomenological formula describes the numerical results quite well, as can be seen 
in figure ITUT b). 



In all the cases considered, we find that the spectral correlations (or more precisely, the 
level repulsion) tend to attenuate fidelity decay. Qualitatively, this can be understood 
from standard perturbation theory JH]; which relates changes in the eigenvector 
components of H to the matrix elements of the perturbation XVap, weighted with the 

o o 

inverse level distance {Ea — Ep)~^. Accordingly, the fidelity decay is slowest in the 
picket-fence case. In the GOE case (which may be associated with chaotic dynamics) it 
is somewhat faster, while in the Poisson case (characteristic of integrable dynamics) it is 
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fastest. In this sense, we can confirm that fidehty is more hkely to be maintained, if the 
system's dynamics is chaotic rather than integrable. This, at first sight counter-intuitive 
result, has been obtained previously by a number of authors [Ulllllllj, although along 
quite different lines of argument. 

3.3. Fidelity distributions 

In particular for eigenstates of Hq, we saw a dramatic effect of the spectral correlations 
on the decay of the average fidelity. As it turned out, sufficiently strong correlations can 
practically stop the decay of {F(t)) at values quite close to one (cf. figure |H1). However, 
the large sample-to-sample fluctuations in this case make us suspect that for individual 
systems this stabilisation effect might be reduced, or not even observable. However, 
in the case of the fidelity amplitude, we could actually prove that if the IPR of the 
initial state is only small enough, then effects of the spectral correlations are observable 
even for individual systems. We expect that the situation will be similar in the case of 
the fidelity itself. There will be a trade off between quantum ergodicity and the size 
of spectral correlations. Besides of its own interest, we study the fidelity distributions 
P{F{t)) to clarify this point. We generate the required histograms at certain equidistant 
time instances. For each histogram, the interval (0, 1] is divided into 50 boxes. Then, 
for a certain sample of rirun systems we count the frequencies with which the fidelity 
for an individual system, at a given time, falls into one of those boxes. After proper 
normalisation the histograms are lined up along the time axis, which results in three 
dimensional plots. 

In figure ^2 we focus on the qualitative features of the fidelity distribution for 
different choices of the initial states. For simplicity we show the Poisson case, though 
note that our observations are equally valid for the GOE case. In all three panels we 
show the fidelity distribution as it evolves in time, together with the behaviour of the 
average fidelity (thick dashed line in the xy-plane). For i^o-eigeiistates, panel (a), the 
maximum of the distribution remains very close to one, while its average continues to 
decay (very slowly, though). In this case, linear response predicts a purely linear decay, 
which can be observed at small times. However, for larger times t > 0.3 the average 
fidelity starts to saturate, an effect we cannot describe theoretically (as we noticed 
already in figure |H1 exponentiation does not help, here). The separation between the 
average and the maximum of the fidelity distribution is due to the development of a 
long range tail. 

For random initial states, panel (b), the fidelity distribution is rather symmetric 
around its average value. Its shape looks almost Gaussian. We checked that the width 
of the distribution decreases roughly proportional to 1/ ^/N. This means that we have 
true quantum ergodicity for the fidelity of a random initial state. For sufficiently high 
dimension of the Hamiltonian, the sample-to-sample fluctuations can be made arbitrarily 
small. 

In panel (c), we show the fidelity distribution, where, for each individual system. 
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Figure 11. Fidelity distribution for N = 100, A = 0.1 for the Poisson case. Histogram 
as a function of time (solid lines), for initial eigenstates (a), for random initial states 
(b), and for the average over all possible initial eigenstates (c). The distributions are 
normalised, such that 6.17. a,t t ~ the height of the bar is 50, because the widths of 
the histogram boxes is 1/50. The average fidelity is plotted as a thick dashed line in 
the xy-plane. 

we have averaged over all available i^o-eigenstates as initial states. This also reduces 
the sample-to-sample fluctuations drastically. Though, qualitatively, the width of the 
fldelity distribution behaves similar to the case with random initial states, it is noticeably 
smaller. Of course, the average fldelity shows here exactly the same behaviour as in 
panel (a). 

In flgure we show examples of fldelity distributions for the Poisson, the GOE, and 
the picket-fence case, at three different times ti, t2, and ^3. At a time ti = 0.195 much 
before the Heisenberg time, at Heisenberg time ^2 = 0.974 and at a time much after the 
Heisenberg time ^3 = 1.884. The different panels are arranged into columns, where in 
the left column we used i^o-eigenstates as initial states, in the middle column we used 
random initial states, and in the right column, we averaged over all i^o-sigenstates. 

In the left column we clearly see the development of a long range tail towards low 
fldelities, which causes the separation between the average and the maximum of the 
fldelity distribution. Though the tail is more pronounced in the Poisson case, it is still 
important in the GOE case, but may be somewhat less important in the picket-fence 
case. In the middle and the right column, the fldelity distributions have almost perfect 
symmetric shapes, at least as long as the distributions do not come too close to the 
boundaries zero and one. In the case where we averaged over the ifg-sigenstates (right 
column) , the widths of the distributions are small enough to allow to distinguish even 
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Figure 12. Fidelity distributions for N = 100, A = 0.1 for the Poisson case (solid line), 
the GOE case (dashed line), and the picket-fence case (dotted line) at three different 
times: t = 0.195 [(a),(b),(c)], t = 0.974 [(d),(e),(f)], and t = 1.884 [(g),(h),(i)]. In the 
left column [(a),(d),(g)], the initial states were eigenstates oi Hq, in the middle column 
[(b), (e), (h)] we used random initial states, and in the right column [(c),(f),(i)] we 
averaged over all initial eigenstates of Hq . The vertical lines give the average value of 
the fidelity for the whole sample for the Poisson case (solid line), for the GOE case 
(dotted line), and for the picket-fence case (dashed line). 



individual systems with respect to the different spectral correlations. For random initial 
states (middle column), the widths of the fidelity distributions are not small enough, 
and one would have to use bigger Hamiltonian matrices (A^ > 1000). 

4. The fidelity in the long time limit 

Here we will calculate the survival probability, the fidelity amplitude, and the fidelity in 
the limit t oo (with the strength of the perturbation held fixed). One may think that 
in general, only the eigenvectors are important. As will be shown below, this is indeed 
the long as Hq is not constructed from a picket-fence spectrum. We use time 

independent perturbation theory, to obtain approximations to the eigenvectors. The 
details can be found in [Appendix B| 

Remember that we are using units for time and energy such that the Heisenberg 
time and the mean level spacing are both equal to one (see section 12)). In these units, 
the time dependent Schrodinger equation reads: 

idt^it) = 271 {Ho + XV) ^{t) = . (30) 

Suppose we can diagonalise H = Hq + XV exactly. Then, an initial state x will evolve 
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in time according to 

u{t) = O A{t) X , A(t) = diag(e-2"'^"*) , H O = diag{E^) . (31) 

We assume that Hq is diagonal, so that we obtain for the evolution in the interaction 
picture: 

x{t) = Ao{ty^ u{t) = Ao{-t) O A{t) X , Ao(t) = diag(e-'"'-^'^*) , (32) 

o 

where Ea are the eigenvalues of Hq. Here, Ao(— t) O A{t) is the exact echo operator 
introduced in equation (^. We expect, that in the long time limit, all dependencies on 
the eigenvalues of H and Hq will drop out, so that the problem consists in finding an 
appropriate approximation for the eigenvectors, i.e. for the orthogonal matrix O. 



The fidelity amplitude can be expressed as the expectation value of the echo operator 
with respect to the state x. 

fit) = (Ao(-t) O Ait) 0^)o = e^-^-^* O^p e'^^^^^* 0,^3 x, . (33) 

In this expression, the phases in the exponential can never cancel. Therefore, the fidelity 
amplitude, averaged over the perturbation V and over Ho, vanishes in the limit t — > oo. 

lini {fit))Qy = . (34) 
For the fidelity, we obtain a more complicated expression: 



Fit) 



e2-^"*0„^ e-^-^/^'O^^ X, 



afd-y 



2 



^ 5^ 5Z 5Z e2'^'(^"-^-')* Oo^fS Oa'(3 O^p Oy^ °X^ Xy . (35) 

aa' 77' 

Here, the arrow indicates that we discarded already those terms containing phases of 
the form exp[— 27ri(i?/3 — Ep/) t], as they cannot contribute to the long time limit of the 
average fidelity. 

If X is an eigenstate of Hq, xp = S^p, then the fidelity as a function of time coincides 
with the survival probability, as we have seen in the beginning of section El For the 
limit value of both quantities, we thus get: 

= lim(F(t))o,y = lim(^(t))o,y = ^^(0^)0,^ = (ipr(0^x))o,y . (36) 

This is the inverse participation ratio of the local density of states (i.e. the projections 
of the i^Q-basis state onto the eigenbasis of H). By contrast, if x is a random state, 
(F(t))o,y ^ F^, where: 



= E E O^, xl)^ ^ + E E e'-'^-^-'' 01, 01, 

Q^7 /3 ' 0^0' /3 ^ 
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a/3 "-^7/3 



= N(WT2) 5 + ' + 3 E O*.) ■ (37) 

For long times, the time dependence survives only in the picket-fence case. Then, at 
integer values of t, the fidelity reaches exactly the value of the survival probability. 
In between, the exponentials give zero on the average. In all other cases, where the 
spectrum of Hq is sufficiently random, the exponentials give always zero. To proceed, 
we will now calculate the limit value of the fidelity (disregarding the exponential). In 
this case: 



oo 



'-^[' + l^(Ot,)oy^ , (38) 



N 



due to X^Q^O^gO^g = 1. Hence, both results are linearly related: 

= ■ (39) 

Perturbation theory for F^ 

Here we consider the long time limit F^ of the fidelity for initial i/o-eigenstates, as 
defined in equation (j36|l . For the GOE well as for the Poisson case, the 

standard (Rayleigh-Schrodinger) perturbation theory (TH] leads to divergences. To 
obtain meaningful results one had to use degenerate perturbation theory, as in |13j . 
Unfortunately, this is quite involved in our case, where powers of matrix elements occur 
up to order four. Therefore we restrict our analytical treatment to the picket-fence 
case. We assume that the initial state is taken from the centre of the spectrum. As the 
spectrum is deterministic, only the average over the GOE perturbation matrix V has 
to be performed. All averages required are of the form (0|^0^^)y. They are calculated 
in [Appendix B[ using the standard perturbation theory. Up to forth order we obtain: 

{0L)v = l-2\'^ + X'^ {OlJv = J^^,^ (40) 
(it is assumed that a (3) and therefore 

i^:> = l-2A^y + 4A^^ + 0(A^). (41) 

Figure ^1 shows F^ as a function of the perturbation strength A for the three 
different spectral ensembles under consideration. In general the (asymptotic) behaviour 
of F^{X) for small A strongly depends on the degree of level repulsion in the respective 
spectrum. Our numerical results range from an apparently linear decrease in the Poisson 
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Figure 13. The average fidelity in the large time limit for initial eigenstates of Hq 
{N = 100). Numerical results for large samples ni-un = 50 000, for the Poisson case 
(0), the GOE (□), and the picket-fence case (O ). Perturbation theory, equation ||1T|| . 
for the picket-fence case (solid line). 



case up to a quadratic decrease in the picket-fence case, while for the GOE case we 
obtain some intermediate behaviour. For large A all curves for F^{X) seem to merge 
into a single one. In the picket-fence case we observe that the perturbative result, 
equation (jHJ, describes the behaviour at small A very well. 

5. Conclusions 

In this paper, we have developed and analysed a random matrix model for echo 
dynamics in a chaotic system. The fidelity has been computed in the linear response 
approximation, and a transition from a linear decay to a quadratic one at times of 
the order of half the Heisenberg time has been found. It is easily verified that in 
the limit of small perturbations one obtains quite generally a Gaussian decay. By 
exponentiating our linear response result, we have obtained a theoretical description 
which reproduces the Gaussian decay in this limit. By construction, this procedure 
provides a fair approximation for the perturbative and the golden rule regime. However, 
by comparing with recent fidelity studies for a deterministic chaotic system [12] , we could 
show that our approximation is accurate in the crossover regime also. 

The main new ingredient in our approach is the dependence of the fidelity decay 
on the spectral two-point form factor of the underlying unperturbed Hamiltonian. This 
lead us to systematically study the effects of spectral correlations on the fidelity decay. 
Besides the GOE type fluctuations which we expect in classically chaotic systems, we 
also investigated random spectra (no correlations) and picket-fence spectra (maximal 
correlations) as the two extreme cases. 

In general, we have found that spectral correlations {i.e. level repulsion) tend to 
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inhibit fidelity decay. However, the effect is usually quite small. It becomes big only 
if we choose eigenstates of the unperturbed Hamiltonian as initial states. Though the 
correlation effects are really dominant in this case, quantum ergodicity is then lost. This 
means that the sample-to-sample fluctuations are as large as the fidelity signal itself, so 
that it is absolutely necessary to perform an ensemble average over many systems, or 
at least many initial states. In order to restore quantum ergodicity, we need to consider 
initial states with large IPR's (if expanded in the unperturbed basis). Unfortunately, 
this inevitably leads to much smaller correlation effects. However, these effects are still 
observable, even for a single fidelity experiment, provided that the IPR is sufficiently 
small (~ 10-3). 

Finally, we considered the fidelity in the limit of large times. We have shown that 
the result for a random initial state is related in a simple way to the result for an 
eigenstate of the unperturbed Hamiltonian. Again we investigated the effect of spectral 
correlations on this limit value for the fidelity. Using perturbation theory, we could 
obtain an analytical result for the picket-fence case, valid for small perturbations. 
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Appendix A. The correlation integral for different statistical spectra 

In order to calculate the correlation integral, equation ()14j) . we need the second integral 
of the two-point form factor. If we denote the first integral by B{t), we may write: 



For a given statistical spectrum {E^} the two-point form factor b2{t) is defined in the 
equations ()1()|1 and (fTTjl . For the GOE and the CUE case, the two-point form factors 
are well known ^3], and for the Poisson case, b2{t) is simply zero. Due to the lack of an 
appropriate reference, we shall calculate 62 (^) for the picket-fence case. 




(A.l) 
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Appendix A. 1. Picket-fence 

In the picket-fence case, we may assume, that Ea = a. With the definition ^(t) = 

^ o o 

^^^j^ exp(27ri(_E'^ — Ea) t), we may write: 

N N 

e(^)=iV-^^g"^g-^ g = e2-S (A.2) 

a=l /3=1 

which allows to obtain ^ (t) in closed form. Namely, using the relation ^"=o ~ 
(1 -g), we get: 



so that 

2-g-g-i l-cos(27rt) sin2(7rt) ^ ' 

Now, we will show that ^(t) — ^ J2nez^(^ ~ ^) limit ^ oo. As ^(t) is periodic 

with period one, it is sufficient to consider t E (— |, |). Let f(t) be an arbitrary smooth 
function, with support in (— |, |). Then: 



sin^TT Nt) /(O) /-^ sin'(7rA ^t) 



AT- 



i J-AT-i sm^(7rt) A y_ 

dt cos2(7rA^t/2) = - /(O) / ds cos2(s) ^ /(O) . (A.5) 

This proves the supposition, so that we can conclude: 

"^Sit-n) = l + 6{t) -b2{t) b2it) = l-J2Kt-n) . (A.6) 

T/ie yirsi and second integral ofb2(t) The first integral of the two-point form factor is: 



B{t) = [ dr 
Jo 



l-J2S{t-n) 
while for the second integral, we obtain: 



t - [t] [t] = max(n < t) , (A.7) 







'drfi(r) = |-[t](t-M±l) . (A.8) 



Appendix A.2. GOE case 

The two-point form factor for the GOE case reads |15j : 

2t - 1 : < t < 1 

1 + t ln(2t- 1) : l<t 
1 - 2t + t ln(2t + 1) + e{t - 1) [2(t -l)-t ln(2t - 1)] , (A.9) 



62 (t) =t ln(2t + l) - 
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where 6{t) is the unit step function. Its first integral gives: 

B{t)= [ dr 62(r) =t-t^+ / dr r ln(2r + 1) 
Jq Jq 



dr 2r - / dr r ln(2r - 1) 



+e{t-i) 
+eit-i) 



(t-l)(t + 2) - 1/4 , , 
{t - if + ^ ^; ^ ln(2t - 1) 



4 2 
while the second integration yields: 

f dr B{r) = 5(^V2 -tV3) + ^ T {t' - 1/4) ln(2r + 1) + 9{t - 1) 



5(t-l)=^ 3(t-l)^ 
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2 7i 



dr (t^ - 1/4) ln(2r - 1) 



The two remaining integrals give: 

Ai(t) = - / dr (t^ - 1/4) ln(2r + 1) = — / ds {s - 2)s In s 
2 Jo 16 Ji 



1 

16 



(z//3 - l)y^ Iny 



y - 1 y - 1 
9 



. y = 2t+l 



1 /■* 1 r'^*^^ 

Mt) = ^ / dr (t' - 1/4) ln(2r - 1) = — ds {s + 2)s In s 



1 

16 



(y/3 + l)2/' Iny- 



16 A 

- 1 _ y"^ - 1 

9 2 



where we have used that: 

r , , 1/^ Iny - 1 

/ ds s m s = : — 

Ji 2 4 



/ ds^ s In s = 
/i 3 



?/=2t-l 

Iny y'^ — 1 



9 



Appendix A. 3. GUE case 

In this case, the two-point form factor is simply: 

1-t : < t < 1 
: l<t 

We give the result for the second integral directly: 

Jo I t/2 - 1/6 : Kt 



hit) 
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(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 



Appendix B. Time independent perturbation theory 



Consider the Hamiltonian H = Hq + A V^. \i Hq has a picket-fence spectrum, we 
may use the standard (Rayleigh-Schrodinger) perturbation theory to compute low order 
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approximations to the eigenvectors of H. A compact derivation of the perturbation series 
can be obtained with the help of the Greens functions for H and Hq (see reference 



G{z) = ^—, Goiz) = ^— ^ Giz) = [l-XGoiz)V]~'Goiz). (B.l) 

Z — 11 Z — -Og 

o 

Assume that Hq is diagonal in the basis used: {a\HQ [3) = 5^^ Ea, and define \ea) as 
the eigenvector of H, such that: 

o 

{ea\H Cfs) = 6ai3 Ea |ea)^|«) E^^E^ (B.2) 
as A ^ 0. Then we may write the projector onto the eigenstate \ea) as follows: 

Pa = \ea) {ea\ = ^^j^zG{z) , (B.3) 

o 

where the integration path is a simple loop enclosing exclusively the pole E^ of Gq{z) 
and the pole E^ of G{z). Then we may obtain the perturbation series as follows: 

p!^^ = leL"^) (eL")| =Y.X' , Tf) = dz [Go{z) Vf G,{z) . (B.4) 

The fi''^ are operators, and the tilde is used to distinguish them from the matrices T^^^ 
introduced below. Note that so far, the projectors Pi"^ are not normalised, so that an 
approximation to the absolute value squared of a matrix element of O reads: 

where = Ylx^aPx^- As a final step, one may expand this expression in a Taylor 
series in A up to any desired order. 

In what follows, we will calculate the matrix elements of O up to forth order in A 
and then average the results over the perturbation V which is assumed to be a GOE 
matrix. To this end, we first consider the operators T^, defined in ()B.4|) : 

Ti°) = ^y"d2 Go{z) = \a) {a\=pa 

f^J-^ = ^ y"dz Go{z) V Go{z) =paVG'^ + G'^ Vp^ , G', = J2 



a — (3 



W = ^ idz Gc^iz) V Gc^(z) V G^iz 
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PaV G'^V G'^ + G'^V po^V G'^ + G'^V G'^V p^ 



t!^^ = Y.^G',Vf Pa[V G'.rK (B.6) 

fc=0 

For the matrices T^^\ these formulae imply that: T*^") = 1, T*^^) = and T'^'^^ has only 
zeros on the diagonal for all A; > 1. 
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Our aim is it to calculate the average of products of the form (0|^0^^)y. These 
are in fact just the one- vector averages of order four of this particular ensemble of 
(approximately) orthogonal matrices. To do so we first expand the squared matrix 
elements of O in terms of the perturbation parameter A. Up to fourth order, we get 

OL = 1 - E A' ^(') + (A^ Si^^f = 1-X' ^(2) - A^ ^(3) - A^ (^(^) - S!^^') (B.7) 

1=2 

01 = E A'^ T^t' (1 - = A^ + A^ + A^ (t^^ - T^'J S^^^) . (B.8) 

k=2 

Now, the desired one-vector averages are easily constructed. To this end we assume 

o 

that the unperturbed energies form the picket-fence spectrum: Ea = a, and we consider 
the limit oo. Thus, it remains to perform the average over the GOE-matrix V. 

Henceforth, we use simple angular brackets without subscript to denote such an average. 
For a 7^ ^ 7^ /i 7^ a, we obtain: 

{OD = 1 - 2A^ - 2A^ (5f ) - 2A^ - (5f ^)) + A^ (5f ^) 

= 1 - 2A^ (5(2)) - A^ (2 {S^f) - 3 (5f ')) (B.9) 

{0}J = A' (T^f) (B.IO) 
{OL 01) = ( (1 - A^ 5(2)) (a2 t(2) + A^ + A^ [Tg) - t(2) sj^^] ) ) 

= ^'{TtJ) + ^'{{T^i)-2{T^'JS^^^)) (B.ll) 

OJJ = A^ (t(2) T(2)) . (B.12) 
There are in total eight different quantities to average: 

(^JV/^TT^IVt-^ (B.15) 



/y{2)y{2)s / >£ \ ('Q^Q^ 

^ \(«-OM«-/i)V (a-OM«-/^)^ ^ ^ 

^ (a-n2 (a - /i)2 +E (a-^)4 
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^ ^-^ \ \ a — f a — II a — u a — 



« — 4 a — n 



a — u a — t, 



111 1 



a — 4 a — f-i a — u 



a — 4^ 



1 1 

+ 



45 9 15 



E 



{a - /i) (a - 



(B.18) 
(B.19) 



TT 
+ ^ 



3 



(B.20) 
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